Superfluidity versus Bloch oscillations in confined atomic gases 
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We study the superfluid properties of (quasi) one-dimensional bosonic atom gases/liquids in traps 
with finite geometries in the presence of strong quantum fluctuations. Driving the condensate with 
a moving defect we find the nucleation rate for phase slips using instanton techniques. While phase 
slips are quenched in a ring resulting in a superfluid response, they proliferate in a tube geometry 
where we find Bloch oscillations in the chemical potential. These Bloch oscillations describe the 
individual tunneling of atoms through the defect and thus are a consequence of particle quantization. 
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Bose-Einstein condensation [0 (BEC) and superfluid- 
ity Q are basic characteristics of Bosonic quantum gases 
and fluids. While the Bose-Einstein condensate (of den- 
sity no) is a thermodynamic quantity characterizing off- 
diagonal-long-range order, the superfluid density n s de- 
scribes the response to a perturbation in the broken phase 
||. In real quantum liquids, such as bulk 4 He, conden- 
sation and superfluidity appear in unison, but in gen- 
eral one may be realized without the other. E.g., non- 
interacting Bose gases in three dimensions form a conden- 
sate without superfluidity (n s — 0), while confining a real 
quantum fluid to two dimensions destroys the condensate 
(no = 0) but preserves superfluidity. In one dimension 
only superfluidity may survive at zero temperature. 

A major recent breakthrough is the realization of a 
Bose-Einstein condensate in weakly interacting atomic 
gases. The confinement within a trap has important 
consequences for the condensate e.g., in a three 

dimensional harmonic trap the well known ideal gas con- 
densation temperature T BE = 3.3h 2 i> 2 / 3 /m is replaced by 
^ho = 0-94 hujhoN 1 / 3 . Quite notable is the appearance 
of a condensate below = huj^N/ \n2N in an ideal 
gas confined to a (quasi) one-dimensional (ID) geometry 
(here, v, N and m are the bulk density, particle number, 
and mass of the bosons, while wfj Q = U]\\w± with us± and 
u)u denoting the transverse and longitudinal trapping fre- 
quencies). The appearance of such a condensate at 
is characterized by a sharp crossover for an ideal ID gas; 
the broadening of the transition due to interparticle in- 
teractions has been discussed by Petrov et ai. ||. 

First attempts to probe the (bulk) superfluid prop- 
erties of condensed atom gases through a moving laser 
beam have been carried out recently 0; the results on 
the critical velocity are in rough agreement with expec- 
tations deriving from a weak coupling analysis based on 
the Gross-Pitaevskii theory An interesting question 
then arises regarding the interplay of superfluidity and 
enhanced thermal/quantum fluctuations due to dimen- 
sional reduction. In this letter we study the superfluid 
response of (quasi) one-dimensional atomic gases and flu- 



ids trapped within finite tube and ring geometries, see 
Fig. 1, where quantum phase slips tend to destroy super- 
fluidity. While interesting on their own, these questions 
have attracted much attention recently through the novel 
atom chip technology (^] allowing for the experimental re- 
alization of strongly confined atom gases exhibiting large 
quantum fluctuations. 




FIG. 1. Trap geometries of the atomic gas perturbed by a 
moving laser beam: (a) ring structure with periodic boundary 
conditions, (b) finite length (L) tube with closed ends. 

The destruction of dissipation free flow in one dimen- 
sional (ID) superconductors and superfluids is triggered 
by the appearance of quantum phase slips as discussed 
by Zaikin et al. fill for metal wires and by Kagan et al. 
|llfl for superfluid rings. Here, we model the action of 
the laser beam through a moving impurity (velocity v) 
and derive a low frequency effective action describing the 
dynamics of the phase difference across. The quantum 
nucleation rate for phase slips determines the response: 
at finite temperature the infinite system exhibits a linear 
response and hence is not superfluid, A/i cx v with A/i 
the drop in the chemical potential across the impurity. 
This contrasts with the ring geometry where interactions 
quench the phase slip nucleation below a critical veloc- 
ity, establishing a superfluid response. In a finite tube 
the quantum phase slips proliferate and the new non- 
superfluid ground state exhibits Bloch oscillations in the 
chemical potential difference across the moving impurity, 
Afi cx shi(2Tmvt) with n the ID atom density. The phys- 
ical origin of these oscillations is found in the particle 
quantization: the moving impurity enhances the particle 
density in front, producing a chemical potential differ- 
ence across the impurity, which in turn is released each 
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time an atom tunnels through the impurity. 

Trapped ID atomic gases have attracted much inter- 
est recently [p|Jl2|-[l^] . Their nature is conveniently de- 
scribed in terms of the dimensionless parameter 7 -1 = 
nh 2 /mg = nl\/2a. Here, n and m denote the (one- 
dimensional) atom density and mass, respectively, and 
l± = y/h/rruJx is the transverse extension of the ground 
state wave function; assuming a contact potential with 
a scattering length a < l± the interaction acquires a 3D 
character and the interaction parameter takes the form 
g — 2h 2 l 2 _/ma, see for corrections. Small and large 
values of 7 correspond to weakly and strongly interact- 
ing bosons: for weakly interacting gases the energy per 
particle is given by the bosonic expression e B (n) = gn/2 
and the Gross-Pitaevskii equation holds. In the strongly 
interacting situation (7 — > 00) the ID Fermion-Boson 
duality becomes manifest and the energy per particle is 
given by the fermionic expression e F (n) = irh 2 n 2 /6m; the 
implications for the density profile have been discussed 
in ]l2] | and an appropriate modification of the Gross- 
Pitaevskii equation has been proposed by Kolomeisky 
et al. JT3fl , It is the latter limit describing impene- 
trable bosons with strong quantum fluctuations we are 
mainly interested in the present work; the requirements 
for the experimental realization of this so called Tonks- 
Girardeau limit have been discussed in Ref. fll^] , A con- 
venient starting point describing the low energy physics 
in both cases is the imaginary time action for the phase 
(f>(x, t) of the bosonic field tp oc exp[i</>] |p^| , 



c Kh 



h.i 



dr / dx 



c s (d x (/)f + — {d T 4>f 
c s 



(1) 



describing sound modes with velocity c s to be cut-off 
at high energies A. For weak interaction, 7 « 1, the 
sound velocity c s = y/ ng/m, the dimensionless param- 
eter K = 7r/ y/j, and the cutoff A = ng derive from 
the Gross-Pitaevskii equation; quantum fluctuations are 
suppressed in this limit. Increasing the interaction, quan- 
tum fluctuations renormalize the sound velocity c s and 
the dimensionless parameter K in the low energy action 
(Q) : strongly interacting bosons with a contact potential 
g5(x) are described by K w 1 + 4/7 while c s K — Trnh/m 
remains unrenormalized (note that K — > 1 in the strong 
coupling or Tonks- Girardeau limit 7 — > 00; parameters 
K < 1 might be realized for bosons with long range in- 
teractions HI). In Eq. © we assume a flat trapping 
potential along the longitudinal direction. A weak longi- 
tudinal trapping potential can be accounted for by space 
dependent parameters K(x) and c s (x); the resulting de- 
formation of the excitation spectrum will not change the 
main results presented below. Equation (Q) produces the 
T = phase correlator ([^(x) - <t){x')] 2 ) ~ In \x - x'\/K, 
hence phase fluctuations destroy the condensate in the 
infinite system. On the other hand, the logarithmic di- 
vergence of the phase correlator is cut off in a finite trap, 



allowing for the definition of a (quasi-) condensate even 
at finite (but low) temperatures as discussed by Petrov 
et al. H (note that the use of (|l]) requires T < A). 

The superfluid properties of the (quasi-) condensate 
can be probed with a moving laser beam |7j which we 
describe by a (strong) impurity potential V ilnp (z,t) — 
#im P <5(£ — v t) with g imp > g, suppressing the particle den- 
sity locally. In the weakly interacting limit we can in- 
tegrate the Gross-Pitaevskii equation over the impurity 
region and derive a Josephson term coupling the left and 
right parts of the atom gas (the 'leads') MJlTf, 



V(<p) = Ej [1 — cos <p{r)] — hnvip, 



(2) 



with Ej = (K 1 'it)ng 2 / 'g imp , to be renormalized in the 
presence of large interactions. The second term drives 
the phase difference tp across the impurity. Small contri- 
butions from higher harmonics do not modify the results 
below which are dominated by large scales; also, on the 
level of (]l]) such terms are irrelevant in the renormaliza- 
tion group sense flq| . The classical (meta-)stable states 
ipj derive from minimizing V(<p) and fall into the inter- 
vals ifj — 2irj G [— 7r, 7r); on a semi-classical level we define 
the associated ground state \ j) of the j-th well. Next, we 
integrate over the phase dynamics (0) in the leads |l7]] 
and obtain the effective action for the phase difference tp 
across the impurity 



drV [<p(t)] 
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Extended leads produce an ohmic kernel Q{lo) ~ 
with a characteristic time scale r c ~ Kh/Ej for the phase 
dynamics; the finite leads in a trap geometry will strongly 
modify the low frequency part of the kernel Q with dra- 
matic consequences for the response. 

Starting from the classical stationary states <fij describ- 
ing a superfluid system we have to account for quantum 
fluctuations introducing transitions between these states 
which potentially destroy the superfluid response. In- 
deed, depending on the low frequency dynamics encoded 
in the behavior of the kernel Qiui — > 0) the relevant 
(semi-classical) instanton solutions [|l9| of (^) will pro- 
vide us with dissipative phase slips and a linear response 
in the infinite system, coherent hopping and Bloch os- 
cillations in the finite tube, or confinement and hence 
superfluid response in the ring. A crucial element enter- 
ing this analysis is the nature of the quantum variable 
tp itself: while tp £ R is an extended variable if differ- 
ent minima tpj are physically distinguishable, tp G [0, 2ir) 
turns into a compact variable if this is not the case. In 
the following we discuss the superfluid response for the 
different geometries in more detail. 

The infinite system is characterized by an ohmic kernel 
Q(lu) — K\uj\; the effective action (||) describes a particle 
in a periodic potential with damping K/2ir. The tun- 
neling between classical minima leads to the excitation 
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of sound modes rendering the states distinguishable; as a 
consequence the phase tp has to be treated as an extended 
variable pcj ] . The system exhibits a quantum phase tran- 
sition at K = 1 pl| separating a non-superfluid ground 
state with a delocalized phase tp at K < 1 from a super- 
fluid ground state with a localized phase at K > f . The 
finite temperature response at K > I is determined by 
the thermally assisted quantum nucleation of phase slips 
p2j: the corresponding action involves a kink-antikink - 
pair separated by the distance r in imaginary time, 
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Using instantons ||19(] the nucleation rate is given by the 
integral T ~ (1/t 2 ) J +lo ° df exp(— S/%). We recover a 
superfluid response at T = with an algebraic rate 
r ~ v 2K ~ 1 , while the response turns linear at finite tem- 
perature, r ~ vT 2K ~ 2 , as thermally activated quantum 
phase slips destroy the phase coherence across the link. 

The ring geometry (see Fig. [|(a)) introduces periodic 
boundary conditions for the phase, (f>(x, r) = <f>{x + L,t), 
with two important consequences: i) the existence of a 
winding number defines an extended quantum variable 
ip, and ii), the sound modes are quantized and exhibit 
a self interaction due to the compactness of the loop, 
modifying the kernel at low frequencies, 

wL f 2Kc s /L, uj<2c s /L, 

Q{uo) = Kuj coth ~ <^ (4) 

2c s [K\ui\, uj>2c s /L. 

The static potential (TiKc s /irL)p 2 /2 describes the kinetic 
energy of the flow in the ring and is easily understood 
when the static solution <p(x) = tp[l/2 + x/L — &(x)] is 
inserted in ([[]) . The additional potential renders the sys- 
tem superfluid p3[ ] : the new minima satisfy the relation 
(Ej /Tin) sinpj = v — v^pj/n, where the first term is the 
usual flow induced by the motion of the impurity, while 
the second term oc vl — Kc s /nL is due to the static po- 
tential. The absolute minimum at pj with |y— 2vi,j\ < vl 
describes a stable superfluid state with a critical velocity 
v c = vl- Indeed, the T = action for a kink-antikink- 
pair exhibits a linear confinement (we assume K > f), 
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and the nucleation of phase slips is quenched for v < Vl- 
At finite temperatures or large drives \v — 2vl]\ > Vl 
incoherent tunneling processes via thermally activated 
quantum nucleation of phase slips describe the equilibra- 
tion of the ring towards its thermal equilibrium as given 
by the appropriate density matrix. At very large drives 
|t> — 2vl] \ 3> vl the system is far from equilibrium and 
the response resembles that of the infinite wire. 

In a finite length tube the flow vanishes at the 
tube ends providing us with the boundary conditions 



d x (j>(—L/2,T) = d x (f>(L/2,T) = 0. The finite size quan- 
tization of the sound modes in the leads introduces a 
smooth low frequency cutoff in the kernel, 

{ K{L/2c s )u 2 , u) < 2c s /L, 

Q(lu) = Kuj tanh ~< (5) 

2c - 1 K\uj\, u > 2c s /L; 

a simple understanding is provided by inserting the dy- 
namic solution 4>(x, t) = ip(r) [1/2 — Q(x)} into (|I|). The 
massive low frequency dynamics renders the tunneling 
between the minima coherent and using instanton tech- 
niques we can determine the hopping amplitude |24| 



W ~ (h/ Tc ) {c s t c /2L) 



K 



(6) 



between semi-classical states In deriving (g) we have 
assumed L > c s t c , see Q for a discussion of short sys- 
tems. We have to distinguish between 'diagonal' tran- 
sitions without excitations of sound modes and 'non- 
diagonal' ones changing the number of sound modes in 
the leads: as long as the time evolution involves only 
diagonal transitions the states \j) are indistinguishable 
and ip is a compact variable p G [0, 2ir). However, in the 
following analysis it is a matter of convenience to choose 
tp extended and compactify only at the end. 

For Ej/(hc s /KL) > 1 the amplitude Ej is larger than 
the plasma frequency of the well and we can study the 
system response within a tight binding analysis. The ac- 
tion (|^) can be transcribed into the Hamiltonian 



where the last term a e = 2i:hnv describes the driving 
force (p = 2nJ2j m the site basis \j)). Applying 

the unitary transformation U = exp [— 2ir i j\f(t)tp] with 
Af(t) — nvt + Ao we eliminate the drive through a redef- 
inition of the hopping amplitudes, 



.}. (8) 



This Hamiltonian is equivalent to that of an electron in 
a crystal driven by an electric force eE = hnv described 
by the vector potential eA/c = eEt = nvt: at zero drive 
the energy eigenstates form a Bloch band, while a finite 
electric field leads to 'Bloch oscillations' (25) (the above 
transformation corresponds to a gauge transformation, 
and the 'quasi-number' Ao accounts for the gauge free- 
dom). At low temperatures and low drives no sound 
waves are excited in the leads and we identify |0) = \j) 
(the compact character of the phase restricts its value to 
a region centered around the potential minimum at ip = 
and tends to establish phase coherence; on the other 
hand, the quantum nature of the phase and the pres- 
ence of phase slips give a finite probability to any value 
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ip E [0, 2ir), thus reducing the phase coherence across the 
impurity, see also j26)). The Hamiltonian (||) reduces 
to H = -W cos[2irAf(t)] |0)(0| and admits the solution 
|0)(t) = exp [-i JdtH(t)] |0). The state of the system 
and its energy depend on the impurity position via the 
'quasi-number' Af] the usual derivatives provide us with 
the chemical potential difference across the impurity 

Afi = {djsfH) = 2nWsm [2nAf(t)] (9) 

and its time evolution due to the drive v 



d t Af = 



(10) 



In the static limit with v — the spectrum maps out a 
Bloch band E(Af) — — W cos(27rAT) (the 'quasi-number' 
Af plays the role of the 'quasi-momentum' k in a pe- 
riodic crystal), while a finite driving force nv leads to 
'Bloch oscillations' in the chemical potential difference 
A/i = 2irW sm(2Trnvt). These oscillations are due to the 
accumulation of particles in front of the impurity, the 
latter allowing only discrete particles to tunnel. Each 
'Umklapp' process describes a particle tunneling through 
the impurity. The behavior of the tube then is dual to 
that of the classic Josephson junction, as is easily seen 
when replacing Af by the phase drop $ and A/i by the 
supercurrent /: the relations (|^) and (10) are equiva- 
lent to Josephson's famous relations I = I c sin <E> and 
3t$ = 2eV/Ti with I c and V the critical current and volt- 
age across the junction. 

For high temperatures and drives processes involving 
frequencies larger than c s /L induce nondiagonal transi- 
tions which compete with the diagonal ones. 'Bloch oscil- 
lations' then disappear above the crossover temperature 
Tl = hc s /L and the critical drive vl/K. At high tem- 
peratures T > Tj, or high drives u>ii£ / K all processes 
are fast and we recover the physics of the infinite wire 
with incoherent tunneling via the quantum nucleation of 
phase slips (we assume K > 1). 

The quantum nucleation of phase slips leads to a trans- 
fer of energy to the bosonic system at high drives v > v l 
but well below the mean field critical velocity Ej/Tin. 
Then the macroscopic quantum tunneling of the phase 
can be observed via the heating of the sample, in analogy 
to the experiment by Raman et al. j7j] (note that our work 
predicts a dissipation free low-drive response and the ap- 
pearance of a critical velocity for both topologies, ring 
and tube). On the other hand, the Bloch oscillations at 
low drives constitute a macroscopic quantum coherence 
phenomenon leading to density fluctuations within the 
leads. Using a second laser beam to probe the oscillating 
densities in the leads allows to measure these fluctua- 
tions, at least in principle. However, as each 'Umklapp' 
process involves only one particle tunneling through the 
impurity these oscillations will be small, thus requiring a 
high sensitivity in the experiment. 

In conclusion, geometric confinement of the atom gas 
boosts the importance of fluctuations. The superfluid 



response strongly depends on the particular geometry: 
in a ring the phase difference across an impurity is well 
defined and the response remains superfluid below the 
critical velocity vl oc 1/L, while in a tube phase slips 
proliferate and driving the system induces 'Bloch oscilla- 
tions' in the chemical potential across the impurity. 
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